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Abstract 



The OS and PDS renormalization schemes for the effective field theory 
with nucleons and pions are investigated. We explain in detail how the renor- 
malization is implemented using local counterterms. Fits to the NN scattering 
data are performed in the 1 So and 3 5i channels for different values of [ir. An 
error analysis indicates that the range of the theory with perturbative pions 
is consistent with 500 MeV. 
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I. INTRODUCTION 



Effective field theory is an important tool for studying nuclear interactions. To describe 
low energy processes involving nucleons and pions in a model independent way all possible 
operators consistent with the symmetries of QCD are included in an effective Lagrangian. 
The short distance strong interaction physics is parameterized by contact interactions, which 
can be thought of as arising from massive states which have been integrated out. The real 
power of effective field theory is that theoretical errors can be estimated in a systematic way. 
Different contributions to an observable are organized by a power counting. This means we 
have a way of organizing the theory as an expansion in Q/A, where Q is a momentum scale 
which characterizes the process under consideration, and A is the range of validity of the 
effective theory. To a given order in Q/A, only a finite number of operators in the effective 
Lagrangian need be retained, and observables can be predicted in terms of a finite number 
of experimentally determined parameters. Theoretical uncertainty in the calculation can be 
reliably estimated and reduced by calculating higher orders in Q/A. Though A is not known 
a priori, it is expected to be set by the masses of the lightest particles which have not been 
explicitly included in the effective Lagrangian. 

In an effective field theory, ultraviolet divergences must be regulated and a renormaliza- 
tion scheme defined. The ultraviolet divergences give a constraint on the power counting 
since when a divergent loop graph occurs one must include a contact operator that can absorb 
the divergence at the same or lower order in Q. This is familiar from pion chiral perturbation 
theory. The choice of regulator cannot affect physical results, but may make implementing 
a renormalization scheme easier. The renormalization scheme and power counting are also 
tied together. In a natural scheme, the renormalized coefficients of the operators in the 
Lagrangian are normal in size based on dimensional analysis with A. Once a power counting 
is established one can translate between different renormalization schemes at a given order 
in Q without changing the physical predictions. 

Counting powers of Q/A in the nuclear effective theory is a subtle issue because of the 
large S-wave scattering length, a. Usually in an effective theory (e.g., chiral perturbation 
theory for pions), the coefficients of operators of dimension n+4 are assumed to scale as A~ n . 
It is then straightforward to examine an arbitrary graph and determine its power in Q/A. 
Applying this approach to the nucleon-nucleon interaction does not work because the large 
scattering lengths introduce an unnatural length scale. Since l/a( 1 S , ) = — 8. 32 Me V and 
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l/a( 3 S'i) = 36.4 MeV, it is necessary to sum corrections that scale like {Qa) n to all orders. 
Therefore, the power counting must be modified. A detailed discussion of power counting in 
the presence of a large scattering length is given in Refs. |1|-§J. In Ref. [|I|], it is pointed out 
that the large scattering length changes the power counting for graphs with intermediate 
nucleons and four-nucleon couplings with no derivatives. In Ref. Kaplan, Savage, and 
Wise (KSW), point out that the effects of the large scattering length can be incorporated 
into the theory by assigning a power counting to the coefficients of four-nucleon operators. 
The power counting for coefficients of operators mediating S-wave transitions is changed, as 
well as other coefficients because of angular momentum mixing. A more detailed discussion 
of this power counting is left to the next section. A similar power counting is discussed in 
Ref. i. 

Two different calculational techniques for the effective theory of nucleons are used in the 
literature. In one approach, the power counting is applied to regulated N-nucleon poten- 
tials and the Schroedinger equation is solved [|l],^|-[7|] . Solving the Schroedinger equation is 
equivalent to including all ladder graphs with the potential as the two-particle irreducible 
kernel (see, for e.g., |J). The second approach, advocated by KSW, is like ordinary chiral 
perturbation theory in that the power counting is applied directly to the Feynman graphs 
which contribute to the amplitude. Here the method for dealing with nuclear bound states 
and infrared divergences that occur at zero kinetic energy is similar to the methods used in 
Non-Relativistic QED and QCD M. A non-relativistic propagator is used which includes the 
kinetic energy term to regulate the infrared divergence. In the Feynman diagram approach, 
dimensional regularization is the most convenient regulator, and analytic results are readily 
obtained. In the potential method, the Schroedinger equation is usually solved numerically. 
In practice, divergences are regulated and renormalized couplings are defined using a finite 



cutoff scheme. In Ref. [fL0|| , it has been explicitly shown that without pions the potential 
method can deal with large scattering lengths, and gives an expansion in Q/A. 

An important aspect of the KSW analysis is the use of a novel renormalization scheme, 
power divergence subtraction (PDS). In PDS, loop integrals in Feynman graphs are regulated 
using dimensional regularization, and poles in both d = 3 and d = 4 are subtracted. The 
subtraction of d — 3 poles gives a power law dependence on the renormalization point, fiR, 
to the coefficients of four-nucleon operators. Let us use C*2 m to denote the coefficient of 
a four-nucleon operator with 2m derivatives (where, for the moment, we restrict ourselves 
to operators mediating S-wave transitions). Choosing ji R ~ Q, graphs with an arbitrary 
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FIG. 1. The leading order contribution to NN scattering in the KSW power counting. 

number of Co(/xr) vertices scale as 1/Q and must be summed to all orders, as shown in 
Fig. |l|. This is precisely the set of graphs that sums corrections that scale like (Qa) n . Higher 
order contributions form a series in Q/A. In Ref. ||, it is emphasized that it is possible 
to phrase the power counting in a scheme independent manner. The choice of scheme 
is simply to give natural sized coefficients which make the power counting manifest, i.e. 
C-zmiHR) ~ 47r/(MA m /^ +1 ), where C2m{HR) are the renormalized couplings. PDS is one 
example of such a scheme. In Ref. ||11]| , it is shown how the KSW power counting can be 
implemented by solving the Schroedinger equation in a finite cutoff scheme. 

Pions can be added to the effective field theory by identifying them as the pseudo- 
Goldstone bosons of the spontaneously broken chiral symmetry of QCD. All operators with 
the correct transformation properties are added to the effective Lagrangian. This includes 
operators with insertions of the light quark mass matrix and derivatives, whose coefficients 
are needed to cancel ultraviolet divergences from loop graphs. In dimensional regularization, 
these ultraviolet divergences are of the form p 2n m 2m /e. For instance, for nucleons in the 3 £i 
channel, the two loop graph with three pions and a two loop graph with two pions and one Co 
have ultraviolet divergences of the form p 2 /e. This pole must be cancelled by a counterterm 
involving a four-nucleon operator with 2 derivatives. Because divergences of the form p 2n /e 
must be cancelled by local counterterms, pion exchange can only be calculated in a model 
independent way if higher derivative contact interactions are included at the same order that 
these divergence occur fll2,l3|. In Weinberg's |l] power counting, pion exchange is included 



in the leading order potential. Therefore, graphs with arbitrary numbers of pions are leading 
order, while the counterterms necessary to cancel the ultraviolet divergences in these graphs 
are subleading. However, the potential method can still be used. As higher order derivative 
operators are added to the potential the accuracy is systematically improved, because the 
onset of the model dependence of the pion summation appears at higher order in Q/A. For 
example, the cutoff dependence of the two pion graph with one Cq will be cancelled by cutoff 
dependence in C 2 . 

Different estimates of the range, A^, of an effective theory of nucleons with perturbative 
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pions exist in the literature. Some authors [|14|-[T6[] argue that A n is as small as m n , so that 
including perturbative pions is superfluous. One estimate of the range is given by KSW 
who conclude that A„- ~ 300 MeV. They point out that in PDS the renormalization group 
equation for the coefficient Co(//r) is modified by the inclusion of pions in such a way that for 
fJ-R ~ 300 MeV, Co(/xr) scales like n R instead of /i^ 1 . Since the power counting is no longer 
manifest above this scale, KSW conclude that the effective theory breaks down at this point. 



In Ref. |L5[ different renormalized couplings are obtained. Here a breakdown of the power 
counting for Ci(jlr) at ~ m n is observed. A crucial question is whether a breakdown 
in the running of the coupling constants is a physical effect or simply an artifact of the 
renormalization scheme. It is dangerous to draw conclusions based on the large momentum 
behavior of the coupling constants because the beta functions of the couplings are scheme 
dependent^. In Ref. ||17|| , a momentum subtraction scheme is introduced where the power 
law dependence of the coupling constants persists even in the presence of pions, and for all 
values oi > 1/a. This scheme is called the OS scheme, since in a relativistic theory it 



might be called an off-shell momentum subtraction scheme. In Ref. \TB\, a similar scheme 
is applied to the spin singlet channel in the theory without pions, where it is shown to give 
results identical to the PDS scheme. The OS scheme is a natural scheme that works with 
arbitrary partial waves and with pions. Thus, the range of the validity of the effective theory 
is not limited by the large fiR behavior of the couplings. PDS is still a useful scheme in which 
to calculate observables. If one splits Cq(hr) into a non-perturbative and perturbative part, 
Co([1r) = Cq(hr) + Cq P ({1r), then Cq P (^r) ~ 1/[1r for all jiR > 1/a. Once this split has 
been performed, it is straightforward to establish relations between the OS and PDS schemes 
order by order in perturbation theory, and any prediction for an observable will be identical 
in the two schemes up to the order in Q/A n to which it is calculated. Since in both schemes 
there is no scale where the power counting breaks down, it is possible that > 300 MeV. 

Physically, one expects the effective theory to be valid up to a threshold where new de- 
grees of freedom can be created on-shell. For elastic nucleon scattering, the relevant physical 



threshold is production of A resonances which occurs at p = J Mn(M& — M^) = 525 MeV 



1 This is in contrast with dimensionless coupling constants like g in QCD. In that case the 
first two coefficients of the beta function are scheme independent, so conclusions based on 
the behavior of the running coupling constant at small coupling (e.g., asymptotic freedom) 
are physical. 
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(the S-wave channels couple only to the A A intermediate state so p = J 2Mn(Ma — M^) = 
740 MeV jn|). Above this scale, the A must be included as an explicit degree of freedom. 
Below this scale, the A can be integrated out leaving an effective theory of pions and nu- 
cleons. Rho exchange becomes relevant at a scale, p ~ m p = 770 MeV. There is also a 
AT* (1440) A" intermediate state with a threshold of p — 685 MeV. One might expect A,,, to 
be of order these thresholds. However, there is an intermediate scale of 300 MeV associ- 
ated with short distance contributions from potential pion exchange^. Using dimensional 
analysis, a graph with the exchange of n + 1 potential pions is suppressed by p/300MeV 
relative to a graph with n potential pions. Comparison of the size of individual graphs is 
scheme dependent (for example the size of graphs differ in MS and in MS). The 300 MeV 
scale applies only to a subset of graphs, and may change once all graphs at a given order in 
Q are included in the estimate. Therefore, 300 MeV can be taken as an order of magnitude 
estimate for the range of the theory, but the actual range may be enhanced or suppressed 
by an additional numerical factor. 

This then motivates the important question: How does one determine the range of the 
effective field theory? This is obviously a question of great practical importance. Theoretical 
arguments can only give an approximate estimate for the range. A good example comes from 
SU(3) chiral perturbation theory. In this strong coupling theory, it is natural to expect that 
the range of the theory is the chiral symmetry breaking scale A x ~ 2^/2 f n = 1200 MeV 
||20|| . However, the convergence of the momentum expansion will depend on the particular 
process under consideration. For instance, in 7r — 7r scattering the range of the expansion 
is set by the threshold for p production, m p = 770 MeV. In this paper, the range of the 
two nucleon effective theory will be estimated using nucleon-nucleon scattering data. Our 
results are consistent with A^ ~ 500 MeV. As we will explain in section VI, the error 
analysis is applied to S rather than to p cot 5 as in Ref . |L5]| . This range does not depend 
on the value of the renormalization point chosen, and is found in both the OS and PDS 
schemes. However, only next-to-leading order calculations have been used so it is hard to 
estimate the error in this value. When higher order corrections are computed, it should be 
possible to obtain a reasonably accurate estimate of the range of the two nucleon effective 
field theory with perturbative pions. This 500 MeV estimate is based solely on the phase 



2 The phrase "potential pion exchange" will be used for a perturbative pion with energy 
independent propagator. This is sometimes called static pion exchange. 



6 



shift data. Predictions for the deuteron electromagnetic form factors are also in reasonable 



agreement with the data ||21|| , but do not probe momenta greater than 200 MeV, and so do 
not provide interesting information regarding the range of the theory. 

In section II, we review the power counting method of KSW 0,0], and the PDS scheme. 
The importance of being able to count factors of the large nucleon mass in a non-relativistic 
effective field theory is discussed. We review the OS scheme, which is compatible with the 
KSW power counting. We describe the procedure for defining the renormalized couplings 
using local counterterms for each of these schemes. 

In section III, we discuss the theory with only nucleons, where A ~ m^. Local coun- 
terterms for both the PDS and OS schemes are computed. These counterterms are used 
to obtain the beta functions for the four-nucleon operators, and we explain why the beta 
functions for the most relevant operators in this theory are one- loop exact. 

The theory with nucleons and pions is analyzed in section IV. In the 3 Si channel, there 
are corrections to the PDS beta functions at all orders in Q. As examples, we compute 
the PDS beta functions for Cq Si \[jlii) to order Q, and for C\ Si \[1r) to order Q°. In this 
channel, even in the limit — > 0, there are logarithmic divergences (poles of the form p 2 /e 
in dimensional regularization) . In the OS scheme, the 3 Si beta functions can be calculated 
exactly. We compute the exact beta functions for C (/iR), C 2 (/z#), and C 4 (/zr) in the OS 
scheme in the ^So and 3 Si channels. In section V, the counterterms for the coupling constant 
D2i.HR) are derived in the OS and PDS schemes. 

In section VI, we discuss why it is important to have fiR independent amplitudes order 
by order in the expansion. In the OS scheme amplitudes are fiR independent, while in PDS 
Hr independent amplitudes can be obtained by treating part of C (/zr) perturbatively. If 



this is not done then the sensitivity to fiR is larger than one might expect [14], for reasons we 
explain. Fits to the data are presented for different values of fiR and the coupling constants 
in both OS and PDS are shown to evolve according to the renormalization group equations. 

In section VII, an error analysis similar to a method due to Lepage || is used to investi- 
gate the range of the effective field theory with perturbative pions at next-to-leading order. 
Weighted fits are performed for the scattering data in both the ^So and 3 Si channels. Our 
results rule out A^ ~ m n , and are consistent with A„- ~ 500 MeV. 
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II. POWER COUNTING AND RENORMALIZATION SCHEMES 



In this section, the KSW power counting and compatible renormalization schemes are 
discussed. The theory containing only nucleon fields is considered first. Some notation 
is set, and we explain why the large nucleon mass does not affect the power counting. 
The renormalized couplings are then defined in terms of local counterterms, and the KSW 
power counting for coefficients of four- nucleon operators is reviewed. Next, we consider the 
theory including pions. We review the power counting for potential and radiation pions, and 
explain the origin of the 300 MeV scale associated with potential pion exchange. The PDS 
renormalization scheme is then discussed and we introduce the OS momentum subtraction 
scheme, which is also compatible with the power counting. 

Below the scale m n , the pion can be integrated out, leaving a theory of non-relativistic 
nucleons interacting via contact interactions. The Lagrangian in the two nucleon sector is 
given by: 

c NN = rt[id t + ^ /(2M) + ...]n-J2 £c2e>£, (2.1) 

s m=0 

where M is the nucleon mass, and the ellipsis refers to relativistic corrections. is an 

operator with 2m spatial derivatives and four-nucleon fields, N. We will work in a basis in 
which the operators mediate transitions between ingoing and outgoing states of definite total 
angular momentum. Our notation is set up to agree with Refs. . The superscript s will 
give the angular momentum quantum numbers of these states in the standard spectroscopic 
notation, 2S+l Lj. If we denote the incoming and outgoing orbital angular momentum by 
L and L f , then any operator mediating a transition between these states must contain at 
least L + L' derivatives. For states with 5 = 0, \L — L'\ = 0, while for states with 5=1, 
\L — L'\ = or 2. For L = L' = the first few terms in the series are 

^ (2-2) 



s.m 



a 



+ 



= C^i^P^NYi^P^N) - ^L- (^P^NYi^P^ V 2 N) + h.c. 

8 l 

(s) 

where the matrices Pj project onto the correct spin and isospin states 

P? So) = -4 (ia 2 ) (iT 2n ) , Pf 3l) = 4= (iWi) iir 2 ) . (2.3) 



The Galilean invariant derivative in Eq. ( |2.2| ) is V ' = V - 2^7 ■ ^ + ^ , and the ellipsis 
denote contributions with more derivatives and/or higher partial waves. 
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The C2m appearing in Eq. (|2.1| ) are bare parameters. To renormalize the theory, the bare 
coupling is separated into a renormalized coupling and counterterms as follows: 

oo 

/^ybare ^finite ruv,M /^yfinitc /~f / , \ \ " rri/~f ( , \ In A\ 

°2m — °2m — b 2m , ^2m ~ ( ~'2m{fJ'R) ~ 6 L '2m\^R) ■ 1^-4) 

n=0 

Note that we divide the counterterms into two classes. The first, which have the superscript 
uv, contain all genuine ultraviolet divergences. These include 1/e poles, if dimensional regu- 
larization is used, or powers and logarithms of the cutoff if a hard cutoff is used. We will also 



include some finite constants (e.g., the — 7 + ln(47r) that is subtracted in MS) if this proves to 
be convenient for keeping expressions compact. By construction, these counterterms are 
independent, but will depend on Cf™ te . The renormalized coupling is denoted C 2m (/ii?). The 
remaining counterterms, S n C2 m {f^R), contain no ultraviolet divergences and will be referred 
to as the finite counterterms. The choice of the finite counterterms differentiates between 
the schemes in our paper. An infinite number of finite counterterms are needed because 
an infinite number of loop graphs are included at leading order. The renormalization is 
carried out order by order in the loop expansion. The superscript n indicates that S a C2 m 
is included at tree level for a graph with n loops. When higher loop graphs are considered, 
the 5 n C2m counterterm takes the place of n loops |22]. For example, at three loops we have 
three loop diagrams with renormalized couplings at the vertices, two loop diagrams with 
a S 1 C counterterm, one loop diagrams with either one 5 2 C or two 5 1 C , s, and a tree level 
diagram with 5 3 C. Examples are given in section III and appendix A. 

For the nucleon theory, the kinematic part of the power counting is very simple [[1IJ21J . Q is 
identified with a typical external momentum characterizing the process under consideration. 
For instance, in elastic nucleon-nucleon scattering Q ~ p, where p is the center of mass 
momentum^. Each nucleon propagator gives a Q~ 2 , each spatial derivative a Q, each time 
derivative a Q 2 , and each loop integration a Q 5 . 

For non-relativistic nucleons, the scale M is much larger than typical momenta. The 
reason for using a non-relativistic expansion is that each graph will scale as a definite power 
of M. To see this, rescale all energies, q° — > q°/M, or equivalently time coordinates, t — > Mt, 



3 For the scattering N(q + p) + N(q — p) — > N(q + p') + N(q — p') it is useful to define 
p = v 'ME - tot — q 2 + ie, where E tot is the total incoming energy, and M is the nucleon 
mass. To simplify the notation we will work in the center of mass frame, q = 0, where 
p 2 = p 2 = p' 2 = ME, and E is the center of mass energy. For external particles, one can 
always translate between E and p using the equations of motion. 
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so that dimensionful quantities have the same size. Since the measure d 4 x ~ M the Lagrange 
density C ~ 1/M. In coordinate space our nucleon fields scale as N(x) ~ M° (in momentum 
space C ~ M and N(p) ~ \fM), so from Eq. (^TJ) C 2m ~ With the M scaling for the 

couplings determined we can find the scaling of any Feynman graph. A nucleon propagator 
gives one power of M, each momentum space loop integration a 1/M, and the M's from 
external lines are cancelled by M's from the states. For graphs that have insertions of the 
four- nucleon operators Np = Nl+Nv — 1, where Np, Nl, Ny are the number of propagators, 
loops and vertices. Thus, at leading order in the non-relativistic expansion any graph built 
out of the interactions in Eq. Q2.2p scales as M _1 since Np — Nl — Ny = — 1. Therefore, 
the 2-^2 scattering amplitude A ~ 1/M. With the definition of A used here this scaling 
gives a finite cross-section in the M — > oo limit which is physically sensible. Since all graphs 
scale the same way with M, M is irrelevant to the power counting. Relativistic corrections 
are included perturbatively [H,^3|], and are generally suppressed by Q 2 /M 2 relative to the 



leading contribution to an observable. This type of correction will not be considered here. 

In the theory with only nucleons, the only graphs relevant to 2 — > 2 scattering are 
bubble chains. Consider a graph Q with L loops in the non-relativistic limit. In dimensional 
regularization, each loop will give a factor Mp/An, and there are L + l vertices, each giving 
a factor — 2Cf" te p 2m . If the operator 2m appears n m times in the graph (L + l = J2m n m) 
the result is: 



g = T7 II A ) V 3 , where j-^2mn m + L. (2.5) 

M m=0 V 47T / m=Q 

If one matches onto the effective range expansion in MS one finds Cf" te ~ Atc a m+1 / (M A m ) 
||. We again see that all graphs Q are proportional to 1/M, which is also true in PDS. The 
large S-wave scattering lengths enhance the importance of some graphs compared to the 
p power counting. This affects the power counting for S-wave couplings, and through the 
mixing, couplings with L and/or V — 2 and S — 1. For other channels we have the usual 
chiral power counting of p's. The power counting for insertions of four-nucleon operators is 
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C^ L,) ^R)P 2m ~ Q q(s ' m) , where 
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m- 1 for L = II = 



g(s, m) 



m for 5 = 1 and (L, L' = 0,2), or (L, V = 2, 



0) 



(2.6) 



= < 



m + 1 for 5 = 1 and L, L' = 2,2 



2m for all other 5", L, and L' 



With the coefficients C 2m scaling as in Eq. (|2.6|), the graph scales as 



where z = ^n m g(s, m) + L . 



(2.7) 
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Note that the power of Q is less than or equal to the power of p : i < j. A useful mnemonic 
for this power counting is 1/a ~ Q, however, the power counting is still valid for Qa 3> 1. 

This Q power counting will be manifest in any renormalization scheme in which the 
C , 2m(/ i i?) scale with n R ~ Q in such a way that Eq. (|2.6| ) is true. At leading order the 
counterterms 5 n C2m(fJ>R) will have the same Q scaling as the coefficient C 2m (/i_R). These 
schemes may differ by contributions in C2 m (n>R) that scale with a larger power of /zr/A, 
since this will not change the power counting at low momentum. 

Let us now discuss the theory with pions. To add pions, we identify them as 
the three pseudo-Goldstone bosons which arise from the breaking of chiral symmetry, 
SU (2) L x SU {2) R — > SU {2) v . With the pions included in this way, we are doing an expansion 
in m 7T /A 7T and p/A^. Note that in this theory, no matter how small p is made the expansion 
parameter will never be smaller than m n /A n . This theory still includes the four-nucleon 
operators in Eq. ( |2.2| ), but the short distance physics parameterized by the coefficients Cim 
is different because the pion is no longer integrated out. In the pion theory, the short dis- 
tance C*2 m coefficients should be independent of the scale m n . All the m n dependence is now 
contained explicitly in powers of the light quark mass matrix in the Lagrangian. 

Pions will be encoded in the representation, E = £ 2 = exp (2iU/f), where 



and / = 130 MeV is the pion decay constant. Under SU(2)l x SU{2)r the fields transform 
as S — > LUffl, £ — > L£U> = U£K>, and N —>■ UN. The chiral covariant derivative is 




(2.8) 



jjfj, _ _|_ ^£Q/i£t _|_ £t0M£)_ With pions we have the following Lagrangian with terms 
involving 0, 1 and 2 nucleons: 
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= 4-Tr (0"E 9 M Et) + ^_ Tr(m g S + m^) 
o 4 

+ ^ jvv,^ - + & (w +^)n (2.9) 

- £ eg Ofi - D^Tcim^P^i^P^N) + . . . . 

Here = |(£m ? £ + £'mq£'), m q = diag(m u , m^) is the quark mass matrix, m 2 = w(m u + 
m<i) where w is a constant, and Qa = 1-25 is the nucleon axial- vector coupling. The ellipsis 
in Eq. (|2.9| ) denote terms with more derivatives and more powers of 7775. 

With pions there are additional complications to the power counting which are 

similar to those encountered in Non-Relativistic QED and QCD ||. The complications arise 
because there are two relevant energy scales for the pions, £% ~ Q 2 /M for potential pions, 
and E-x ~ Q for radiation pions. When the energy integral in loops is performed via contour 
integration, the graphs with potential pions come from terms in which one keeps the residue 
of a nucleon propagator pole. In these loops, the energy of the loop momentum is ~ Q 2 /M 
and the energy dependent pieces of the pion propagator are suppressed by an additional 
Q 2 /M 2 . Nucleon propagators give a Q~ 2 and the loop integrals give Q 5 . There are also 
radiation pion graphs, in which the residue of the pion pole is kept. In these loops, the loop 
energy is ~ Q, the nucleon propagators give a Q -1 , and loop integrals give a Q 4 . We find 
that adding a radiation pion to a bubble chain of contact interactions gives an additional 
suppression factor Q , where N is the number of nucleon propagators in the radiation pion 
loop. In either case, a pion propagator gives a Q~ 2 , and each nNN vertex gives a Q. The 
combined propagator and vertices for a single pion exchange give Q°, so the pions can be 
treated perturbatively 0. 

In general, pion exchange gives both long and short distance contributions. The short 
distance contributions from potential pions are important since they may limit the range of 
the effective field theory. A single potential pion exchange gives 

■ 9a g • gag Q ■ gyg = -9a_ 

If 2 q 2 + ml 1 2 2f 2 

where the spin indices connect to nucleon fields N^NpN^Ng which belong to external lines 
or propagators. The first term dominates for q 2 ^> m 2 , and can be isolated by taking the 
limit m w — > 0. Graphs with radiation pions are suppressed by m w /M. In the non-relativistic 
limit, with only potential pions, the only loop diagrams are ladders. Consider an arbitrary 
graph Q with n m four point vertices, Cf"* 13 , and k potential pions. For L loops, this graph 



q-a af3 q- a 1& m%q-a aP q- a l5 



Q ' 



Q ' 



+ mi 



(2.10) 
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has a total of L + 1 = k + J2 m n m vertices, and with m„ = 



«/ m=0 J m=0 

oo 

where j= ^] 2 m n m + L . (2-11) 

m=0 

In the ^So channel, the relation in Eq. ( |2.11| ) becomes an equality. The graph Q ~ Q 1 where 
i is given in Eq. ( |2.7| ). The power counting of the 5 nv C2 m counterterms is determined by the 
need to cancel ultraviolet divergences, and will not spoil the scaling for the renormalized 
coefficients, since i < j. For graphs with only potential pions (n m = 0), it appears that our 
expansion is in pj (300 MeV) since 

0f -(SOOMeVr 1 . (2.12) 

Comparing the size of potential pion graphs therefore predicts a range of 300 MeV, but the 
size of these graphs may change depending on the renormalization scheme (i.e., the finite 
subtractions). It is not known a priori how the contact interactions will affect the range of 
the effective theory. The scale 300 MeV is therefore an approximate estimate for the range 
of the effective field theory with perturbative pions. A further discussion of this issue will 
be taken up in section IV. 

Next, consider the power counting for coefficients that multiply operators with powers of 
m q . If we are interested in momenta of order m n , then one counts m q ~ ~ Q 2 . Therefore, 
any interaction term that has an operator with a total of 2m powers of p and will scale 
as 

Qq{s,m) where q ^ 

m) is given in Eq. (|2~lf). For example, D2 m\ ~ Q°. It is important 
to understand that in the KSW power counting D 2 should be treated perturbatively even 
though the structure of the operator it multiplies is similar to that of the leading 4 nucleon 
operator with no derivatives. Graphs with radiation pions will also give contributions with 
powers of m\. In Ref. the power counting for a general radiation pion graph is worked 
out; the only change is the power counting for the coefficients in Eq. (|2.6| ). 



A. PDS 

PDS is one scheme in which the KSW power counting is manifest. In PDS, we first let 
d = 4 and take the S nv C2 m counterterms to subtract 1/e poles as in MS. We use the notation 
for the renormalization point, and \x for the dimensional regularization parameter. In 
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PDS, like in the MS scheme, one takes p, = pr. In a momentum subtraction scheme this 
is not necessary. The next step in PDS is to take d = 3 and define the finite counterterms, 
S n C2m([iR), to subtract the l/(d — 3) poles in the amplitude. Graphs which contribute are 
those whose vertices have a total of 2m derivatives. When calculating the S n C2 m {pR) we 
take m n = since, for instance, counterterms proportional to m\ renormalize coefficients like 
£*2 (/-£/?)• After making these subtractions everything is continued back to four dimensions. 
It is this second set of finite subtractions that gives the right power law dependence on /jr. 
To define the coefficients that multiply operators with powers of m q , a similar procedure 
is followed except we count the powers of at the vertices. In PDS with just nucleons, 
all the graphs that affect the running of C2 m (/iR) are order Q q ^ s ' m \ except for those with 
intermediate states of different orbital angular momentum. For example, the beta function 
for C[ Sl ^ has contributions ~ Q (g( 3 Si,4) = 1), as well as contributions ~ Q 3 from graphs 

rig _3£)i ) 

with two C*2 vertices. When pions are included there are additional graphs that are 

sub- leading in the power counting and effect the running of the couplings. In fact, in section 
IV we will show that there will be corrections to the PDS beta function for Cq Si \pr) at all 
orders in Q. 



B. OS 

Another renormalization scheme that can be used to reproduce the power counting in 
Eq. Q2.6p is a momentum subtraction scheme. A simple physical definition for the renor- 
malized couplings can be made by relating the couplings to the amplitude evaluated at the 
unphysical momentum p = ipR. This scheme will be called the OS scheme, since in a rela- 
tivistic field theory this would be referred to as an off-shell momentum subtraction scheme. 
We start by dividing up the full amplitude as 

oo 

iA' = i E^ + --- ■ (2-13) 

m=0 

Here A\ m contains the Feynman diagrams that will be used to define the coupling C^mi^R) 
(or equivalently the counterterms S n C2 m ). The ellipsis in Eq. (|2.13j ) denotes pieces that 
vanish as m n — > which are not needed to define Ci m (/Ar)- Aj m is defined to contain the 
remaining graphs that scale as Q q ^ s,m \ where q(s,m) is defined in Eq. ( |2.6| ). The definition 
for the renormalized coupling is then 
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P = l[L R "\ \P = l\l R 



- X + XX + XXX + ••• + co — 

FIG. 2. Renormalization condition for Cq^h). The amplitude iA s, ° is the full four point 
function with Cq S \hr) and 8 n CQ ((j,r) vertices, evaluated between incoming and outgoing states 
s = 1 S or 3 5i. 

iA s 2m \ . =-4lMfe) 2m . (2-14) 
v = 

= 

As we will see, this ensures that C2m{nn) scales in the desired way. In general, there may 
be divergent graphs scaling as Q l and p 2m [i < 2m) whose 1/e poles need to be absorbed by 
a 5 uv C*2 m counterterm. For example, consider the graph with two pions and one Cq shown 
in row four of Fig. |5]. This graph has a p 2 /e pole which is cancelled by a counterterm 6 UV C 2 - 
The finite part of this graph is used in Eq. ( |2.14| ) to define C 4 (/^) because the graph is order 
Q. The key point is that since q(s,m) < 2m, an ultraviolet divergence that appears in a 
graph of a given order can always be absorbed into a coefficient that appeared at the same 
or lower order in the power counting. Therefore, we will define 5 uv C 2m in MS to subtract all 
four dimensional 1/e poles so that these subtractions are independent of the renormalization 
point. The finite counterterms are then fixed by the renormalization condition in Eq. ( f2.14|) . 



In the OS scheme, the coupling Co({Ir) is defined by the renormalization condition in 
Fig. H] (with or without pions). This condition is to be imposed order by order in the loop 
expansion so that graphs with n loops determine <5 n Co(/ii?). The m n = part of pion graphs 
will also contribute to C*2 m (/iR) for m > 1 in which case the condition m^ = in Eq. (p. 14 ) 



is important. In the theory with pions, we also need to define couplings multiplying powers 
of m q , like D 2 in Eq. ( |2.9| ). To define these couplings we will not include all the terms in 
the amplitude proportional to m\. In particular, pion exchange graphs give long distance 
non-analytic contributions which will not be used to define the running of the short distance 
coupling D 2 (fiR). The idea that long distance physics must be excluded from the short 
distance coefficients is discussed in Ref. fT5| . A detailed discussion of how we define D 2 (hr) 
in the OS scheme will be left to section V. 
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Note that in the OS scheme there is another approach for calculating an amplitude in 
terms of renormalized couplings. One can calculate all loop graphs in ^4| m in terms of the 
finite (or MS) parameters and then demand that the renormalization condition in Eq. (|2.14|) 
is satisfied. This gives expressions for the renormalized couplings in terms of the constants 
Cf° it0 . The amplitude can then be written in terms of renormalized couplings by inverting 
these equations. This simplifies higher order calculations. 

In the OS scheme, when an amplitude is written in terms of renormalized couplings it will 
be explicitly (jlr independent at each order in Q. The /ir dependence in PDS with pions is 
cancelled by higher order terms. It is possible to obtain hr independent amplitudes in PDS 
if part of Co(fiR) is treated perturbatively jLTj]. Consequences of this fiR dependence will be 
discussed in section VI. In section III, we will see that for the theory with just nucleons the 
OS scheme gives very similar definitions for the renormalized couplings to those in PDS. In 
section IV, we investigate the running couplings in both schemes in the theory with pions. 



III. THEORY WITH PIONS INTEGRATED OUT 

In this section, we compute the renormalized couplings in the non-relativistic nucleon 
effective theory without pions. We expect A ~ m n . This theory will be examined in 
both PDS and the OS scheme. The renormalization program is implemented by explicitly 
calculating the local counterterms. In Ref. [I8|, it is shown that the PDS and OS schemes 



give the same renormalized coupling constants in the 1 S'o channel. Here we also consider 
the spin-triplet channel and higher derivative operators. Divergences in loop integrals are 
regulated using dimensional regularization. For the OS scheme, the same renormalization 
program can be carried out using a momentum cutoff regulator as shown in appendix A. 
Following Ref. [[|, we will multiply each loop integral by (fx/2)( 4 ~ d \ and define d = 4 — 2e. 
Since there are no logarithmic divergences in the nucleon theory, 5 uv C*2m = in dimensional 
regularization. 

In both PDS and OS scheme, it is straightforward to derive the finite counterterms, 
<5 n C 2m (/i_R). The tree level graphs with Cq(hr) and C 2 (/iR) satisfy the renormalization 
condition in Eq. fl2.14p . Therefore, in both PDS and OS, 5° Co = S°C2 = 0. At one and two 
loops we have the graphs in Fig. |3|. In d dimensions, the two graphs in the first row give 

HCo)2 (^)r (¥) (|) 4 -Vz^)- +i ^ , (3 .D 
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2>CX X 

^^2 ^ ^^2 ^ 2 ^ 2 

^jxx X 

FIG. 3. One and two loop counterterms for Co and C*2. The solid lines are nucleon propagators, 
and symmetry factors are shown explicitly. The generalization to higher loops is straightforward. 



determining 5 1 Cq. In PDS, we define the counterterm to cancel the d = 3 pole in Eq. ( p . 1|) 
and then continue back to four dimensions. In the OS scheme, we take d = 4 and demand 
that the contribution to the amplitude in Eq. (|3.1| ) satisfies the condition in Fig. |2|. The 
counterterms calculated in each scheme are the same (with [i = fi R in PDS). In both schemes 
the counterterms determined from the graphs in Fig. |3] are 



^CoW = (^) C M 2 , = -(^) CoGM 3 , (3.2) 

5 l C 2 {^ R ) = 2(^) CaMCbW - ^O^b) = -^j'^NCoW 2 . 



Note that it is essential that loop graphs also have vertices with insertions of the countert- 
erms. For instance, the contribution to the amplitude from all the graphs in the second row 
of Fig. |3| is 

_ !Co (^p^)\ (3.3) 

If the one-loop graph with a 5 1 Cq counterterm had been left out then the answer would 
have been proportional to (p 2 + fi 2 R ) which is not correct. Since the loops in the nucleon 
theory factorize, the renormalized n-loop graph gives {ip+/J* R ) n - Loop graphs will not always 
factorize once pions are included. 

It is straightforward to extend this calculation to n loops and to include higher deriva- 
tives. In both the OS and PDS schemes, this gives the following counterterms (s = ^So, 3 Si, 
n > 1): 
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(3.4) 



n+l 



*-cf i} w = (-ir 1 (^) cf 

TOf^M = (^)" (n + l) ^W<#*ViO, 

n+l ^^V^'M-^-^N 



V 47T 
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.n+l 



47T 
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Note that with fiR ~ Q, the counterterms have the same Q scaling as their corresponding 
coupling constant. In the PDS scheme, there are also subleading terms that come from the 
mixing of angular momentum states. In PDS 



5 n ct Sl \m) 



l fMn R \ n ( 3 Sl) 



47T 



.( 3 Si)/ 



(3.5) 



where the last term is suppressed by Q 2 . In the OS scheme 



S n cf Sl \fM R 



-1 



,n+l 



4tt 



(3.6) 



which is the same as the Sq channel. In the OS scheme, graphs with two C\ couplings 

(3g \ (3g \ 

and any number of Cq s contribute to the beta function for Cg since they are order 
Q 3 . One might also ask about channels where the large scattering length does not effect 
the power counting. In this case C^mi^R) ~ Q°, an d we recover the usual chiral power 
counting. In our OS scheme, the counterterms 8 n c!£l(ii R ) in these channels are either zero 
or a constant independent of fi R . 

From Eq. ( |2.4j ) one can derive the beta functions using 



d 00 d 

film = ^RT, C*2m(/ i i?) = / , l^RT\ ^CWt^i?) 

dm 9fM R 



(3.7) 



The first few beta functions are 



\2 



3 Si, 3 Bi: (3.8) 

in agreement with Refs. [0,0- For S* = states the beta functions are one loop exact in the 
sense that the contribution in Eq. fl3.8| ) comes from the one-loop graphs, with the higher 
order graphs giving contributions which cancel. The reason for this cancellation is that the 
only loop corrections are in the bubble chain, and they form a geometric series. The sum 
of bubble graphs is just the chain of irreducible one loop bubbles for the full (point-like) 
propagator. An analogy would be QED, if the only possible graphs were the two point 
photon graphs with electron loops. In this case the beta function would also be one-loop 
exact because the graphs that are not 1PI do not contribute. In general, the beta functions 
of higher order couplings may have contributions beyond one-loop in cases where angular 
momentum mixing is present. 

Expressions for the running coupling constants can be derived by summing the coun- 
terterms in Eq. ( |2.4| ) or by solving renormalization group equations. For s = 1 Sq or 3 5*i this 
gives 



Cq (hr) - — ? jj7— , C 2 S (n 



qhmte ^ 



Q finite 

M/ir-i ' v - ! >' ' ' " " ^finite \ 2 



i _ Mm. 



2 ' 



(3.9) 



where C£ nite and Cf nite are constants which can be determined by specifying boundary 
conditions. Since the theory should be good for arbitrarily small momenta, one possibility is 
to demand that the amplitude reproduces the effective range expansion, pcot (5) = — 1/a + 
\tqP 2 + 0(p 4 ). In Refs. [0,0] this matching was done at [Mr — giving Co mte = 
(jftmte _ fam. o£a^ ^y e cou ](j equally well have chosen a different matching point, and 
obtained the same results. For ijlr ~ Q, the running couplings in Eq. ( |3.9| ) have the scaling 
in Eq. ( f2.6| ). Written in terms of renormalized couplings the amplitude in the ^So or 3 Si 
channels is B 
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(3.10) 



and satisfies Eq. ([2.14J) . The amplitude A is hr independent. It is interesting to note that 
we can choose a renormalization point where all loop corrections vanish giving 



A s 



J2A 

4:71 
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2m 



m=A 

An 

~Ml/a + ip ~ M 



—ip) p 



2 m 



1 



1 



r 2 , 
- P +.... 



(3.11) 



1/a + ip j 

The amplitude exactly reproduces the effective range expansion by construction. From 
Eq. (|3.11|) the range of the effective field theory can be estimated as A ~ 2/r ~ m n as 
expected. 

It is possible to choose the boundary condition for Co(^r) to change the location of the 



pole that appears at each order in the expansion. For processes involving the deuteron ||2T|J24 
a more natural boundary condition is to choose the pole to appear at — ip = 7 



so Cq 111 * = 47t/ (M7). To recover the effective range expansion, Cq(h) is divided into non 
perturbative and perturbative parts |17| , Cq(hr) = Cq P (hr)+C p (hr), where Cq P (hr) 
and C p ([ir) ~ Q°. In this case the amplitude becomes 



aon- 

1/0 



A s --- 

M 



1 
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47T C^ R ) 



+ 



4vr C 2 (hr) 



P 



(3.12) 



_ 1 + ip M(CZ P (n R ))*( 1 + ip) 2 M(C "%^))*(7 + ^) 
where the first term is order 1/Q, and the second and third terms are order Q°. The RGE's 



are 







d 



Mfi! 
47T 



(3.13) 



^^-C^R^M + OiQ). 



L dfi R uyr ' lw ~ An 

These can be derived by substituting Cq(hr) = Cq P (hr) + C p ([j,r) into the renormaliza- 
tion group equation for Co(/xr). They can also be derived using the counterterm method 
described above. If we demand that the observed scattering length and effective range are 
reproduced at this order then we find 



4vr C p ((x R ) 



1 

7 - - 

a 



An C 2 {hr) 



To 

2 



M(CHa*k)) 2 «' M {Cq P {hr)) 2 

In order for the power counting of C p (hr) to be consistent we must treat 7 — 1/a ~ Q 2 . 



(3.14) 
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FIG. 4. Zero, one, and two-loop graphs with Cq and <5 n Co vertices and potential pion exchange. 
The dashed lines denote potential pion propagators. 



In this section, we study the renormalization of contact interactions in the effective field 
theory with pions. In the 3 Si channel, graphs with two or more consecutive potential pions do 
not factorize and give poles of the form p n /e. We explicitly compute these poles for two loop 
pion graphs. There are also rrv^/e poles in both the 1 £o and 3 Si channels at order Q° [fj||. 
Because of these 1/e poles, pions cannot be summed to all orders in a model independent 
way. The finite counterterms in PDS and OS are different in this theory. Throughout this 
section we will take = 0, since we are only interested in the couplings C 2m (/ii?). The 
-D2 counterterms will be considered in section V. We compute the PDS counterterms 
and beta functions for Co(hr) and C^/ij?) to order Q. In PDS, Co(/xr) no longer obeys 
the Q scaling for fiR > 300 MeV This can be fixed by treating part of the coupling 
C (fiR) perturbatively as discussed in Section VI. The exact expressions for C (//r), C 2 (//r), 
and C&(fiR) are given in the OS scheme and exhibit the correct Q scaling for all fiR > I /a. 
Therefore, it is no longer apparent that the power counting breaks down at 300 MeV. The 
300 MeV scale does appear in the short distance contribution to the amplitude from pion 
exchange, however, it can only be taken as an estimate for the range of the effective field 
theory once pion and contact interactions are both included. In section VI, we will discuss 
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FIG. 5. The basic order Q graphs in the 3 Si channel whose loop integrals do not factorize even 
for m-Tr = 0. 

how experimental data suggests that ^ 300 MeV. 

To determine how the pions contribute to the beta functions for C 2m (//R), we use the 

rules in section II. Some of the pion graphs that will be needed are shown in Fig. |j. 

In both PDS and OS, the first step is to subtract 1/e poles. For two nucleons in the 
channel the spinor indices in Eq. (|2.10|) are dotted into 5 a s <5g 7 . Therefore the = 



piece of pion exchange reduces to a contact interaction and gives no 1/e poles. In the 3 Si 
channel, graphs with two or more consecutive pions do not factorize and may have 1/e poles. 
Order Q graphs with two consecutive potential pions are shown in the first row of Fig. |5|, 
and labeled a), b), and c). We find 

I_2 7 + — -41n(2) + 2h/ ^ 



2 \2pJ V 4vr 

( 9a V f~ipM\ 2 



e 3 V — p 2 — is 

1 -3 7 -61n2 + ^ + 31n^ ^ 



e 6 V — p 2 — ie 



Graphs b) and c) have been written with Cq 111 * 6 vertices to emphasize that the uv counterterm 
which cancels their divergent part is independent of /xr. The divergence in b) is cancelled 
by a tree level graph with the counterterm 



2, uv/i a /nffinite 



S 2 ' UV C 2 = -6C 



2 \ 2 

(4.2) 



— - 7 + In (tt) + 2 - 2 In (2) 



V87T/V 

where the superscript 2 indicates that the counterterm comes in at two loops. The extra 
factor 2 — 2 In (2) is included because this leads to simpler analytic expressions. Expanding 
the Cq bubble graph (second row, first column of Fig. f|) in e gives 
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When graphs with 1/e poles are dressed with Co bubbles, the factors of [2 — 7 — 2 In 2 + In (71)] 
that appear are cancelled by similar factors from the counterterms. In fact, 5 2,nv C2 is the 
only uv counterterm we need for two potential pion exchange with = 0. The 1/e pole in 
c) is nonanalytic since it is proportional to p 3 . When graph c) is added to graphs d) and e) 
the poles cancel. These cancellations continue to occur when more Cq bubbles are added to 
b) and c). After including graphs with 5 2,UV C2 we find 
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(4.4) 
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-p z 
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Note that for fi ~ p there are no large numerical factors from these graphs. 

In the 3 S\ channel, potential pion graphs without contact interactions also have p 2 /e 
poles. The two loop graph with three potential pions (fourth row, second column in Fig. |J) 
is equal to 



V 



3 

- + 

e 



(4.5) 



M \8nf 2 J 

In the Q power counting, this graph is order Q 2 and will not be considered here. Because 
of these 1/e poles it is not possible to sum pion ladder graphs to all orders. Now that the 
ultraviolet divergences have been removed from graphs b) and c), the finite subtractions can 
be performed. 



A. PDS 



For PDS in the x Sq channel, we can compute the effect of potential pions on the C2m{^R) 
counterterms to all orders in Q (neglecting relativistic corrections). For Co(/xr), the relevant 
zero, one, and two loop graphs are shown in Fig. |j. The Co([Ir) and C^/-^) counterterms 
are 
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The PDS counterterms in the 3 Si channel will only be computed to order Q since the 
loop graphs with consecutive pions do not factorize. For this case it is essential to use the 
counterterms to carry out the PDS renormalization program. To define Cq{^r) at order 
Q, we set up the finite subtractions as in Fig. f|, but leave out all graphs with more than 
two potential pions since they are 0(Q 2 ) (we also neglect relativistic corrections that are 
order Q but come with an additional 1/M 2 ). Note that in d = 3 only the overall divergence 
(oc l/{d — 3) n for n loops) is needed since loops with counterterms will cancel the sub- 
divergences. Evaluating the graphs in Fig. || with d = 3 and then continuing back to d — 4 
gives 
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for n > 2 



Note that for graphs with two consecutive potential pions, the /i# dependence does not come 
in the linear combination hr + ip. For instance, adding the PDS counterterm to graph a) 
in Fig. H gives the linear combination 3ip/2 + 9/iR. 

In PDS, like in MS, the renormalized coupling ^(/xr) will depend on \n(fi 2 R /fil) in such 
a way that the ln(/jfj) dependence in the amplitudes in Eq. ( |4.1| ) is cancelled. Here /xo is an 
arbitrary scale expected to be of order A„-. At order Q we find 
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C 2 (a*r) 

for n > 2 . 



Note that the part of <5 n C 2 (/ii?) proportional to In ((jl r / has a coefficient that sums up to 
at this order. From Eqs. (|4.6|) , (|4.8|) , and (|4.9|) we find 
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Note that in the channel all contributions to the beta functions beyond one-loop cancel, 
leaving them one-loop exact. In Ref. ||, the last two terms in (3$ are absent, but should 
be included in the complete order Q calculation. Dimensional analysis implies that the 3 jSi 
beta functions can have corrections at all higher orders in Q, since there is nothing to prevent 
the dimensionless factor (fj,Rg\M)/(87if 2 ) ~ Q from appearing. In Ref. Hl5| , expressions 
for the beta functions are derived by demanding that dA/dfiR = 0, but these are not the 
PDS beta functions. Since in all renormalization schemes dA/dfiR = 0, this condition is not 
sufficient to fix the renormalization scheme uniquely. As discussed in Ref. JT7J, the large hr 
behavior of Cq Si \^r) is unknown because of the higher order corrections. 



B. OS 

In the OS scheme, there is no such ambiguity since at a given order in Q the running of 
all the coupling constants that enter at that order are known exactly. The coupling {/^r) 
has contributions only from the nucleon graphs discussed in section II and therefore has the 
same beta function. For C^^/ir), the order Q° graphs in A2 include the nucleon graphs 
from section II, as well as the graphs with one potential pion and any number of Co vertices. 
At tree level we add a finite counterterm to cancel the m, = part of the tree level pion 
interaction at p = i\ir 

^ s W = -f|4r. (4-11) 

Z J H'R 

This counterterm is order Q~ 2 like C2(//r) itself. Since all the graphs in A2 factorize the 
higher loop counterterms are the same as in the theory without pions, so S n C2 for n > 1 are 
given in Eq. (|3.4j) . The exact beta function is then 

tf> = 2 ^« Ca{m) CM + 2 Ml + ^Cofe)) 2 -L , (4.12) 
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Note that the finite ln(/^ 2 /(— p 2 — ie)) terms in Eq. ( f4.1|) are order Q and in the OS scheme 
do not affect the running of C2(//r), but rather C^r). In terms of the finite constants 



Co mte and Cf nite we have solutions 



/~< finite 
°0 



Mm ' 

(4tt) 



finite 



Cg nitc M 



n 2 • 



Att 



(4.13) 



Although it may seem that the piece of C^^r) that goes as will spoil the power 

counting for low momentum, in fact, the l//^ part dominates entirely until ~ 1/a, since 
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next-to- leading order OS amplitude is 
~C 2 {fi R ) p 2 
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(4.14) 



which is order Q° as desired. 

One might still ask if the problem with the 300 MeV scale will reappear in higher order 
coefficients. To check that this is not the case we compute the running of the coupling 
Ca{hr) in the OS scheme. The easiest way to compute this running coupling constant is to 
compute the order Q amplitude in terms of the finite couplings, Cf" te , and then demand 
that the amplitude satisfies the renormalization condition in Eq. ( |2.14| ). The graphs we need 
to compute include those with 



n 
a) 

Hi) 
iv) 



one C4 and any number of C 's , 

two CVs and any number of Cq's , 

one C2, one potential pion and any number of CVs. 

two potential pions and any number of Cq's . 



(4.15) 



Computing these graphs in terms of the finite couplings and then demanding that the 
amplitudes satisfy the renormalization condition gives the OS couplings 
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where here \x is an unknown scale expected to be of order A„-. Again the pion contributions do 
not spoil the lir scaling behavior, since they are suppressed by factors of the large scattering 
length. Note that at order Q the PDS coupling C^lir) is the Qa — ► limit of Eq. ( |4.16| ). 

In this section, expressions for the renormalized couplings Co(lir), C 2 (lir), and C±(lir) 
were derived in the PDS and OS schemes working to order Q. For the 3 Si channel, we have 
shown that C${lir) has corrections at all orders in Q in PDS. Unlike PDS, the OS couplings 
C 2m {^R) can be computed exactly because they only have contributions at one order in Q. 
The OS couplings exhibit the correct llr scaling for all lir> 1/a. 



V. THE COUPLING D 2 (hr) 

In this section, the OS and PDS counterterms for D 2 (hr) are computed. To define 
D 2 (^r) in the OS scheme, we take 

iA s {D 2 ) . =-zD { 2 s \ Ll R)ml, (5.1) 

P=t/4Rj 

where A S (D 2 ) contains terms in the amplitude that are analytic in m 2 and proportional 
to m 2 . Only terms that are analytic in m 2 are kept because it is unnatural to put long- 
distance nonanalytic contributions that come from pion exchange into the definition of the 
short distance coupling [T5|. For example, one potential pion exchange gives a m^/p 2 ln(l + 



4p 2 /m 2 ) term. Including this in A S (D 2 ) would give D 2 (hr) both a branch cut at lir = m n /2 
as well as explicit dependence on the scale m^. In the OS scheme, D 2 {^r) will be calculated 
as follows. First m 2 / e poles are subtracted. The finite counterterms are then determined by 
including graphs with a single D 2 {^,r) or potential pion and any number of Cq(llr) vertices 
in A S (D 2 ). Contributions from these graphs that are non-analytic in m 2 are dropped. 

There is a m 2 /e pole in the O(Q ) graph in the third row and third column of Fig. |] 
OfJH, so we have a counterterm 



/ M c hmte \ 2 n 2 



j- ~ 7 + log (?r) 



V 4tt ) Af 2 [2e 1 ' ~ bV " y J ' (5 ' 2) 
Note that when this counterterm is dressed with Co bubbles the extra factors of 2 — In 2 
from Eq. ( |4.^| ) will cancel without the need for an additional finite term in 5 2,nv D 2 . After 
subtracting this counterterm the value of the two-loop graph is 
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For PDS we set /i = and then find finite counterterms 
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(5.4) 
(5.5) 
for n > 2 . 



Here /xq is an unknown scale expected to be of order A n . 

In the OS scheme, the 5 1 D2(fi R ) counterterm is the same as in PDS. In dimensional 
regularization logarithms of the form ln(/i 2 /m 2 ) will appear in loop graphs. To make the /i 2 
dependent part analytic in m 2 we write 
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and then subtract the ln(fi 2 /fi R ) term with the counterterms. This will give D 2 (/Xr) a /i^ 
dependence which cancels the ln(/i|,/m 2 ) in the amplitude. In the OS scheme, the m\j2 in 
Eq. ( |5.3|) gets subtracted along with the logarithm. We find 
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Summing the counterterms the solutions for D2(fi R ) are then 
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which can be written as 
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where /i 2 = /1q exp 



/-64vr 2 / 2 D 2 lnit0 
VM 2 ^(C finite ) 2 



/i 2 = /i 2 exp I 1 



647T 2 / 2 D 2 lnitc 



in PDS , 
in OS. 



(5.8) 

(5.9) 
(5.10) 



M 2 ^(C finitc ) 2 

The parameter jl must be determined by fitting to data. With m n ~ Q ~ //r, D2{[i R )m 2 K 
Q° in both OS and PDS, implying that /^(/^ii) should be treated perturbatively. 
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VI. SCHEMES AND AMPLITUDES 



In this section, the amplitudes in the 1 Sq and 3 Si channels are presented to order Q°, 
both in PDS [@,[| and OS. Fits to the 1 S and 3 Si phase shift data are done in both 
schemes for different values of /ir. As pointed out in Ref. |17| , one has the freedom to split 
C ({j,r) into perturbative and nonperturbative pieces: C ([j,r) = Cq P (hr) + Cq(/Ar), where 
Cq P {[Ir) ~ Q~ l and Cq(hr) ~ Q°. This division is necessary in PDS in order to obtain 
(jLr independent amplitudes at each order. Furthermore, Cq P (hr) ~ 1/ Hr so the coefficients 
scale in a manner consistent with the power counting for all fi R > 1/a. For convenience we 
will drop the superscript np in what follows. Some issues that arise in matching the pion 
theory onto the effective range expansion are also discussed. 

First, we give the nucleon-nucleon scattering amplitudes in the PDS and OS schemes. 
In PDS, the amplitudes were calculated to order Q° in Refs. 0,0]. At this order, amplitudes 
in the and 3 Si channels have the same functional form, 

A = + A^ + A^ + OiQ 1 ) . (6.1) 

In both OS and PDS we have 
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The remaining part of the order Q° PDS amplitude is 
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Note that since we have made a different finite subtraction than KSW the second term has 
a prefactor of 1/2, rather than a 1 as in Ref. 0]. In the OS scheme 
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In appendix B, we give relations between the OS and PDS couplings that appear at this 
order. Using these equations it is easy to verify that the amplitudes are equivalent in the 
two schemes. 
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In the PDS scheme, there are order Q° contributions to flo (c.f., Eq. (|4.10|) ). If the order 
Q° contributions are separated from the order 1/Q pieces, the beta function for CqGur) is 

dC^ R ) Mfj, R , . r ^ 



^fl— a = 2 — ; — C (hr) 



C P oM + ff- 2 +0(Q). (6.5) 



This equation has the solution 



«J (6 6) 

where K is a constant which must be order Q 2 for Cq([1r) ~ Q°. (Recall, from Eq. (|3.14 ) 
that K = '-f — 1/a < 1/ain the pure nucleon theory.) Including Cq(^lr) makes the PDS 
amplitudes explicitly fiR independent. In performing fits to the data the constant K and 
coupling D 2 (hr) cannot be determined independently. In what follows we will drop K and 
simply remember that the values of D 2 {[Ir) extracted from the fits may differ from the 
renormalized coupling in the Lagrangian. In PDS, if Cq(^ir) is omitted from our expressions 
then D 2 (hr) does not follow the renormalization group equation, as we will see below. 

In the OS scheme, the constant g 2 A /{2f 2 ) in Eqs.( p.5|) and (|6.6|) is not present, so 
Cq(hr)/Cq(iir) 2 is fiR independent. The OS scheme amplitudes A* -1 ) and A^ ' are therefore 
fiR independent without Cq([ir) as can be seen by examining Eqs. ( |4.13|) , and (|5.7|). In OS 
the constant K will also be absorbed into D 2 (hr). 



Using the Nijmegen phase shifts |25[ between 7 and 100 MeV, we fit the coefficients 
Co(^r), C 2 ([1r) and D 2 (fj, R ). We took m n = 137 MeV. Clearly we would like to bias the fit 
towards the low momentum points since that is where the theoretical error is smallest. This 
can be accomplished by assigning a percent error, ~ p/ (300 MeV), to the data and then 
minimizing the x 2 function. In Tables | and |T| we show the values^ of Co(/ar), C 2 (iir) and 
D 2 ([ir) extracted from the fits for fiR = 70, 100, 137, 160, 280 MeV. These values exhibit the 
[Mr dependence predicted by the RGE's to ~ 1% in the 1 S channel and ~ 4% in the 3 Si 
channel. In Fig. the results of the fits are shown. The results of the fits shown in the 
figure are identical in both schemes. Higher order corrections will give contributions to 5 of 
the form p 2 /A 2 . The error in 5 at p = 300 MeV is consistent with > 500 MeV. 



4 The coefficients extracted from our fits differ from those in Ref. because we have 
emphasized the low energy data as opposed to doing a global fit. It is interesting to note 
that using our PDS value C 2 (fiR = 137 MeV) = 11.5 fm 4 , the prediction for the RMS charge 
radius of the deuteron plj becomes 1.966 fm which is within 1% of the experimental result. 
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Fit to 






Fit to 3 Si 




Mfl(MeV) 




c 2 {m) 






C 2 (iir) 


D 2 {m) 


70 


-6.48 


10.11 


-0.532 


-22.73 


171. 


-70.41 


100 


—4.71 


5.36 


1.763 


-9.93 


32.7 


—4.157 


137 


-3.53 


3.01 


2.000 


-5.88 


11.5 


1.500 


160 


-3.05 


2.25 


1.869 


-4.69 


7.32 


1.897 


280 


-1.79 


0.772 


1.105 


-2.19 




1.UU1 


TABLE I. 


and 3 iSi 


couplings in 


the PDS scheme. Co(/j,r) 


(in fm ), C 2 {fi 


r) (in tm J, and 


D 2 (/j,r) (in fm 4 ) are fit to the Nijmegen 


data at different values of 










Fit to 






Fit to 3 5i 




Mii(MeV) 


Cq(hr) 


c 2 (m) 


D 2 {m) 


Co(hr) 


C 2 (hr) 


D 2 (hr) 


70 


-6.50 


9.75 


-6.047 


-24.1 


121. 


-170.1 


100 


-4.73 


5.33 


-1.143 


-10.0 


27.3 


-20.18 


137 


-3.54 


3.00 


0.378 


-5.92 


10.5 


-4.124 


160 


-3.06 


2.25 


0.658 


-4.74 


6.89 


-1.671 


280 


-1.80 


0.779 


0.692 


-2.23 


1.61 


0.2985 



TABLE II. 1 Sq and 3 S\ couplings in the OS scheme. Cq(^r) (in fm 2 ), C 2 ([ir) (in fm 4 ), and 
D 2 (^r) (in fm 4 ) are fit to the Nijmegen data at different values of hr. 



For processes involving the deuteron it is convenient to fix Cq([jlr) using the deuteron 
binding energy, Co(m^) = — 5.708 fm 2 . With this constraint we find C^m^) = 10.80 fm 4 and 
D 2 (m w ) = 1.075 fm 4 in the PDS scheme. The fit to the phase shift data with these values is 
as good as that in Fig. |6|. 

In PDS, it is necessary to break Cq(hr) into perturbative and non-perturbative parts to 
obtain amplitudes that are fiR independent order- by-order. If Cq^r) is omitted then the 
values of D 2 {[Ir) determined from the fit will not follow the RGE. To see this we define the 
Hr independent quantity 

(6.7) 



R 



-D 2 (jjlr) | MfMg 2 ^ ln //4 



(4 



C q (jjl r )* SttVStt/ 2 , 

and choose the constant c so that R = 1 for [ir = 137 MeV. For other values of [ir the 
deviation of R from 1 gives the discrepancy between the values predicted by the RGE and 
those extracted from the fit. For /j, R = 70, 280 MeV we find R = —0.53, 7.25 in the l S 
channel and R = —0.52, 11.4 in the 3 Si channel. These large deviations disappear if Cq(ijlr) 
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FIG. 6. Fit to the phase shift data emphasizing the low momentum region. The solid line is 
the Nijmegen fit to the data [25|, the long dashed line is the order 1/Q result, and the short dashed 
line is the order Q° result. 

is included. Without CqO^r), the PDS amplitude is still independent to the order that 
one is working. However, as explained below, this residual /ir dependence gives larger 
corrections than expected |14| since it makes the tuning that was setup to give the large 
scattering length hr dependent. 

For momenta p <C the pion can be integrated out leaving the effective range expansion 



P COt (5) = — - + y p 2 + V 2 P 4 + V 3 p 6 + f 4 p 8 + . . . . 



(6.8) 



Performing a matching calculation between the two theories gives expressions for a, r$ and 
the Vi in terms of the parameters in the pion theory. Since the theory with pions is an 
expansion in Q these predictions take the form of Taylor series in Q/A n 

i = 7 + ?.«J\ r ° 

a 



7 + y = ££f\ v n = £ B». ,0.0) 



8=2 * i=0 i=2-2n 

where B W ~ Q 1 . At this time only the first coefficient in each series is known since p cot 5 
has only been calculated to order Q 2 . The notation 
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will be used to denote the location of the perturbative pole in the amplitudes. In PDS 
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0.01468 
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32.57 


0.01101 


137 


-10.16 


2.04 


-8.121 


0.01480 


48.96 


-16.76 


32.19 


0.01098 


280 


-10.23 


2.12 


-8.105 


0.01484 


46.39 


-12.64 


33.76 


0.01111 



(2) 

TABLE III. Values of 7, Bq , l/a, and tq (in MeV) obtained from our fits. Three values of 
Hr are shown to emphasize that the value of the extracted parameters depends weakly on /j,r. 
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Note that if Cq(^r) had been neglected then B$ would not be hr independent. With 
Hr = m n Eq. ( |6.11| ) agrees with Ref. || if their definition of D 2 (fj, R ) is adopted. In the OS 
scheme we have 
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The value of the remaining determined at this order are the same in both schemes 
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For n = 2,3, 4, Eq. ( |6.13| ) gives the low-energy theorems derived in Ref. |]TB[ if we set 7 = l/a. 

Recall that the unnaturally large scattering length a is a fine tuning that was accounted 
for by demanding that in Eq. (|6.10|) , Co(^r) is close to its ultraviolet fixed point, and 7 ~ 
l/a. Examining the expression for l/a in Eq. (|6.9| ) it may seem that this could be destroyed 
by chiral corrections. If D 2 ([ir) ~ Cq(^r) 2 then the first term gives B {2) ~ 205 MeV. In fact 
from Table |T|, we see that the fit gives B {2) < l/a. The reason for this small value is that 
since oc (A^ 1 ') 2 the amplitude has a double pole. Since this pole is spurious (occurring 
from the perturbative expansion) the residue of the double pole must be small in order to 
fit the data. This leads to a good fit condition which will be approximately satisfied 
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As explained in Ref. [17], this condition implies ~ An^ 2 /M. In fact this gives the right 

(2) 

order of magnitude for the values of Bq determined from the fits in Table |H11 Similar 
good fit conditions occur at higher order keeping the coefficients Bq small. Thus the 
tuning 7 « 1/a is not destroyed, but instead naturally kept by the form of the perturbative 
expansion. The division of Cq(hr) into nonperturbative and perturbative pieces is arbitrary, 
allowing us to set up the theory so that the Q expansion for 1/a is well behaved. 

In Table |lJl| we see that when Bq is added to 7, values of 1/a are obtained which 
are close to the physical values, l/a( 1 S , ) = — 8. 32 Me V and l/a( 3 Si) = 36.4 MeV. It is 
encouraging that the value of 7 found from fits in the 3 Si channel are close to the physical 
pole in the amplitude which corresponds to the deuteron, 7 = 45.7 MeV. Values for tq can 
also be predicted from the fits using Eq. (|6.12| ). Experimentally, ro( 1 5o) = 0.0139 MeV -1 
and r ( 3 Si) = 0.00888 MeV" 1 , so the values in Table agree to the expected accuracy. It 
is not yet clear |17| whether values of the Vi extracted from experimental data ||16],|26| are 
accurate enough to test the low-energy theorems for v 2 , v 3 , and V4. 



VII. DETERMINING THE RANGE A,, 

Here we will examine the phase shift data to see what it tells us about the range of the 



effective field theory with perturbative pions. In Ref. [15], a Lepage analysis is performed on 
the observable pcot 5(p) in the 1 Sq channel. Near 350 MeV the experimental 1 Sq phase shift 
passes through zero. Therefore, the error |pcot5 NPWA — pcot5 EFT | is greatly exaggerated 
since p cot 8{p) — ■> 00. To avoid this problem we will use the l So and 3 Si phase shifts as our 
observables, since A5 = |5 NPWA - <5 EFT | remains finite for all p. The next-to-leading order 
amplitudes given in section V will be used. The phase shifts have an expansion of the form 
5 = + 5(1) + 0(Q7 2 ), where § 



5(°) = --ln 
2 



i +i £U<-» 

Itx 



SW = P Ji Ai0) (7 1) 

ZTT 



Recall that a momentum expansion of 5 would result in terms with only odd powers of 
p. However, the expansion for 5 in Eq. (|7.1|) is not simply a momentum expansion, so the 
next-to-leading order calculation can have errors which scale as p 2 /A 2 r . For example, once 
pions are included we can have a term p 2 tan" 1 (2p/m 7r ) which is odd in p, order Q 2 , and 
scales as p 2 for large momenta. 
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FIG. 7. Error analysis for the phase shifts in the and 3 <Si channels. AS is the difference 
between the the effective field theory prediction and the Nijmegen partial wave analysis p5[| . The 
long and short dashed lines use the O(Q ) and 0{Q) theoretical phase shifts respectively. 

In Fig. |7|, we plot AS versus p using log-log axes. Note that the sharp dips in Fig. [7| are 
just locations where the theory happens to agree with the data exactly. The Nijmegen data 



25fl is available up to p = 405 MeV. In a theory with just a momentum expansion the errors 
will appear as straight lines on the log-log plot as pointed out by Lepage ||. In the pion 
theory the expansion is in both m^/A^ and p/A n , so this is no longer true. For p > m n we 
expect the errors to be of the formQ 
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The fact that there is always a p/A n error arises from the fact that, as seen in Eq. ( |6.9|) , r 
is reproduced in the effective field theory as an expansion in m^/A^. For p/A n 3> m^/A^ 
the slope of the lines on the plot should indicate the lowest power of p that has not been 
included. At low momentum the error in AS^ is dominated by the pm"/A™ +1 term and 
the lines should be parallel. From Fig. |7| we see that the error is smallest at low momentum 
and increases as the momentum increases, which is how the theoretical error is expected to 



5 At momenta 1/a C p <C m n we could have AS^ ~ B^/p ~ m^/A^p. However, as 
explained in section VI, B^ < 1/a so this term is very small. 
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behave. 

It is clear that even for p ~ 400 MeV the next-to-leading order calculations are reducing 
the error in the phase shift. Because two new parameters are added at next-to- leading order 
it is always possible to force exact agreement at some value of p. However, if one were to 
force the data to agree too well at high momentum then this would destroy the agreement 
at low momentum. Since the improvement of the fit in Fig. [I] at high momentum does not 
come at the expense of the fit at low momentum this is evidence that the error is being 
reduced in a systematic way. At high momentum one expects that the error is ~ p 2 / 'A 2 . 
From Fig. [7], AS ~ 0.26 radians for p = 400 MeV, implying ~ 800 MeV. This is only a 
rough estimate for the range because we cannot yet exclude the possibility that the next- 
to- next-to-leading order phase shift has an anomalously small coefficient. Even though the 
lines in Fig. [7] are not straight they should still cross at approximately the range of the 
theory since at this point higher order corrections do not improve the agreement with the 
data. This error analysis is consistent with the possibility that the range is > 500 MeV. 

Further information on the range of the effective field theory can be obtained by ex- 
amining electromagnetic processes involving the deuteron pl|j2^| , such as the deuteron 
charge radius, electromagnetic form factors, deuteron polarizability, and deuteron Comp- 
ton scattering. For these observables errors are typically ~ 30 — 40% at leading order and 
~ 10% at next-to-leading order. This is what one would expect if the expansion parameter 
m^/A^ ~ 1/3, implying A„- ~ 410 MeV. This is consistent with our previous estimate for 
A„.. If the range is this large one should expect that the error in deuteron properties will be 
at the few percent level once next-to-next-to-leading order calculations are performed. 

VIII. CONCLUSION 

In this paper the structure of the effective field theories of nucleons with and with- 
out pions is studied. We discuss a momentum subtraction scheme, the OS scheme, which 
obeys the KSW power counting. The method of local counterterms is used to obtain the 
renormalization group equations for the coupling constants in these theories. Using local 
counterterms defines the OS and PDS renormalization schemes unambiguously. Two-loop 
graphs with potential pions in the 3 Si channel are computed and shown to have p 2 /e poles. 
The presence of 1/e poles implies that the only model independent piece of pion exchange is 
the part that can be treated perturbatively. We obtain the renormalized couplings Co(/xr), 
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C*2 and C^/xr) at order Q in the OS and PDS schemes. 

We have emphasized why it is important to have llr independent amplitudes order by 
order in Q. Such amplitudes are obtained automatically in the OS scheme. In PDS llr 
independent amplitudes may be obtained by treating part of Co(lir) perturbatively. Another 
result concerns the large lir behavior of the couplings in this theory. In the OS scheme the 
coupling constants obey the KSW power counting for all llr > 1/a. In PDS the breakdown 
in the power counting for Co(lir) is avoided if Co(lir) is split into non-perturbative and 
perturbative parts. Therefore, the breakdown of the scaling in the coupling constants is 
artificial. 

Next-to-leading order calculations of nucleon-nucleon phase shift data || provide fits to 
data at large momenta which are far more accurate than one would expect if the theory 
broke down completely at 300 MeV. Of course, this does not mean that nucleon effective 
theory can be applied at arbitrarily high energies. The scale, Mg\/(Siif 2 ) ~ 300 MeV, is 
associated with short distance contributions from pion exchange and provides an order of 
magnitude estimate for the range. In the S-wave channel, A production and p exchange 
become relevant at ~ 700 MeV, which sets an upper limit on the range of the expansion. 
To get a better understanding of the range of the nucleon effective theory with perturbative 
pions one must examine experimental data. An error analysis of the S wave phase shifts with 
next-to-leading-order calculations seems to be consistent with a range of 500 MeV. Though 
next-to-next-leading order corrections need to be compared with data and other processes 
investigated, we remain cautiously optimistic that the range could be as large as 500 MeV. 
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APPENDIX A: LOOP INTEGRALS WITH A MOMENTUM CUTOFF 

REGULATOR 



Although the analysis in section III used dimensional regularization to regulate divergent 
loop integrals, the results for the coefficients C2 m (fiR) in our momentum subtraction scheme 
are independent of this choice. As an exercise we will derive the counterterms for Co(hr) 
and CziuFt) using a momentum cutoff regulator, A. This will give us the chance to see what 
type of complications can arise using a different regulator. Note that this is not the same 
as using a finite cutoff scheme. There the momentum cutoff plays a double role as both a 
regulator and as part of the subtraction scheme. 

The graph in the first row first column of Fig. ^ gives 



and the same finite counterterm, 8 1 Co(hr) in Eq. ([3.4Q is used to satisfy the condition in 
Fig. 0. The renormalized graph is then the same as calculated in dimensional regularization 
in section III. Note that contributions of order p 2 have been neglected in defining Cq (jir) 
as required by our renormalization condition. An added complication with a cutoff is that 
graphs with only Co's give a contribution to the amplitude proportional to p 2 . However, as 
A — > oo, p tanh (p/A) — > 0, so these terms can be completely neglected. This will remain 
true even for higher loops since the counterterms will always cancel dangerous powers of A 
that appear in the numerator. At n loops we find an ultraviolet counterterm of the form 




(Al) 



An ultraviolet counterterm cancels the linear divergence, 




(A2) 




(A3) 



while the finite counterterms are given by Eq. ( |3.4| ). 

The graph in the third row first column of Fig. [5] gives 




n «M^ ^ [A 3 2 r mp ,-\ V 

2 — ^C C 2 I hp AH - ptanh 1 - 

2vr 2 \ 6 y [ 2 y \A 



) 



(A4) 
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Note that there are different contributions from this graph when the vertices are in the order 
Co Co C2 or Co C2 Co- At order p 2 , this graph gives a correction to the counterterm 5 1,uv Co, 
i.e., 8 x ' uv Co -> 5^ W C + 5 u > uv C 0l where 

M A 3 

rl*,uv/^y o /nrflmte /"(finite / a r\ 

5 C = -— 2C C 2 — . (A5) 

Unlike the contribution to 5 1,nv Co in Eq. ( |A3| ) , 5 1 * ,nv Co is to be treated perturbatively, so that 
at it only appears once in any graph. The justification of this fact is that this contribution 
to the counterterm appeared at order Q° (a purely formal trick to recover this counting is 
to take A ~ fiR ~ Q). The counterterm 5 l,uv C2 is fixed by considering the order p 2 terms in 
Fig. H row 3. From Eq. (|A4j) (the tanh -1 piece can again be thrown away) we have 

S 1 ' UY C 2 = — 2 C finitc C 2 finitc (-—) ■ (A6) 

4"7T V 7T / 

The calculation for higher loops is similar and there are again corrections 5 n *' uv Co to 5 n ' uv Co 

/ — M\ n / ?A\ n_1 A 3 

5 n *' uv C = n(n + 1) (C finite ) n C 2 finite — , 

\ An J V 7r / 3n 

S n,u VC2 = _(_^_\ ( n + l)( C <finitc ) n C ,finitc (_^M . (A7) 

The finite counterterms are the same as in Eq. ( |3.4| ). Thus the running couplings and 
amplitudes with a cutoff are the same as found using dimensional regularization. 



APPENDIX B: RELATIONS BETWEEN COUPLINGS IN OS AND PDS 

Here we give explicit relations between the coupling constants that occur at order Q° in 
the OS and PDS schemes. Couplings on the left are in PDS while those on the right are in 
the OS scheme. 



Co{^r) 
C 2 (hr) 
C0M 2 

D 2 (H R ) MfMgl 

Co^r) 2 



Co(flR 



m + 9 a 



7T V87T/ 2 
1 



Co(hr) 

c 2 (m) 

Co(^r) : 



+ 



9a 



2 1 



R) 



C (/ij 



R 



V 2 C (m) 



Co(hr) 



9\ r 



2f 



C Q (/ij 



+ 



R 



47T 



Co (A* 



+ 



M/i 



R 



R 



4tt 



(Bl) 



As in section VI the superscript np on Co(hr) has been dropped. 
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